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ABSTRACT 



Context. Clumping in the winds of massive stars may significantly reduce empirical mass-loss rates, and which in turn may have a 
large impact on our understanding of massive star evolution. 

Aims. Here, we investigate wind-clumping through the linear polarization induced by light scattering off the clumps. 
Methods. Through the use of an analytic wind clumping model, we predict the time evolution of the linear polarimetry over a large 
parameter space. We concentrate on the Luminous Blue Variables, which display the greatest amount of polarimetric variability and 
for which we recently conducted a spectropolarimetric survey. 

Results. Our model results indicate that the observed level of polarimetric variability can be reproduced for two regimes of parameter 
space: one of a small number of massive, optically-thick clumps; and one of a very large number of low-mass clumps. 
Conclusions. Although a systematic time-resolved monitoring campaign is required to distinguish between the two scenarios, we 
currently favour the latter, given the short timescale of the observed polarization variability. As the polarization is predicted to scale 
linearly with mass-loss rate, we anticipate that all hot stars with very large mass-loss rates should display polarimetric variability. This 
is consistent with recent findings that intrinsic polarization is more common in stars with strong Ha emission. 

Key words. Stars: mass-loss - Stars: winds, outflows - Polarization 



1. Introduction 

Many different species of massive star dis play intrinsic linear 
polar ization - from main-sequence O stars (Lu pie & Nord sieck 
I1987L LN87 hereafter), to post main-sequen ce transitional 
objects such as OB superg iants (lHavesI U984al LN87), Red 
Supergiants ( Haveill984bl) . Luminous Blue Variables (LBVs , 
iDavies et alJ I2005I) . B[e] sup ergiants dMelgareio et all |2001|) 
and Wolf-Rayet stars (WRs, ISt.-Louisetal.lll987l) . This in- 
trinsic polarization results from the scattering of starlight by 
aspherical structure in the stellar wind. In the case of the 
B[e]SGs, this inhomogeneity is thought to be in the form 
of a circumstellar disk, although the data for this is limited 
(e.g. [Zickgraf & Schulte-Ladbeck|[l9"89t ISchulte-Ladbeck et all 
1993). The intrinsic polarization of the LBV AG Carina was 
also t hought to result from a glob al axi-symmetric wind struc- 
ture (Sch ulte-Ladbeck et al.l 11994b . However, multi-epoch ob- 
servations have shown th at the polarization angle of AG Car 
is stochastically variable (IDavies et al.ll2.Q05b . The polarization 
of P Cygni was found to be variable on timescales of ^ weeks 
dNordsiecketal.1 120011) . which is similar to the short timescale 
variability observed in WRs dRobert et al J l989). This behaviour 
is commonly explained as being due to wind-clumping (e.g. 
LN87). 

Not just a quirk of massive stars, wind-clumping could have 
serious implications for our understanding of stellar evolution. 
Mass-loss rates have traditionally been determined by fitting re- 
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combination lines with models that assume a smooth ou tflow 
(for LBVs see e.g. IStahl et afl 120011: iMachado et al.ll2002h . As 
recombination is a two-body interaction process it depends on 
p 2 , so inhomogeneity will result in an over-estimation of the den- 
sity (as (p 2 )/(p) 2 ), and hence an over-estimation of the mass-loss 
rate. 

Recently, there has been much discussion about clumping 
in the winds of O stars. The value of the clumping factor 
/ci = (p 2 ) /(p) 2 is as yet an ope n is sue. On the on e hand , Ha 
analyses bv lRepolust et all (120041) and lMokiem et all d2007l) find 
that a modest amount of wind clumping, with a clumping fac- 
tor up to / C ] a 5 is required to m atch the smooth line-driven 
wind mo dels of IVink et al.l (120001) . But on the other hand, X- 
ray (e.g.|O skinova et al. 200 6]) and quan titative UV studies (e.g. 
Boure t et alj2005t iFullerton et al.l 2006) call for f c \ values much 
larger than 10. If confirmed, this could significantly affect our 
understanding of massive star evolution. 

In addition to spectral diagnostic techniques, linear po- 
larimetry studies have proven to be a powerful tool to probe the 
geometry of the wind structure. For i nstance, the results of a 
spectropolarimetric survey of LBVs by IDavies et al.l d2005l) call 
for an onset of wind clumping in close proximity to the stellar 
photosphere, as to be able to reproduce the observed levels of po- 
larization. LBVs may be an ideal testbed to study hot star wind 
clumping as the strength of the intrinsic polarization is much 
larger than that observed in O and WR stars. Systematic monitor- 
ing campaigns have shown that WRs display polarization around 
the 0.1% level, with variability on very short timescales of < 1 
day. Additionally, the variability was found to have a greater 
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amplitude for those stars with slower terminal w ind speeds ( 
Voo = 1200 3600kms^, AP ~ 0.1 -> 0.02%; iRobert et al.l 
1 19891) . LBVs have wind speeds around a factor of 10 slower 
and show intrinsic polarization of ^0.5%. Although there have 
yet to be any syste matic polar imetric monitoring campaigns of 
LBVs, the data of Haves] JT985I) appeared to indicate variability 
on timescales comparable to the wind flow-time (ffl = R^/v,*,), 
which for LBVs is on the order of ~days. This is consistent with 
the polari zation arising within a few stellar radii of the central 
star (c.f. ICassinelli etal]|1987b . 

Theoretical studies of clumping-induced polarization are 
few, and tend to be ' first-step' investigations rather than 
comprehensive s t udies. iRod rigues & Magalhaes ( 2000j) and 
ICode & Whitnevl d 19951) presented Mon te-Carlo studie s of the 
polarization produced by a single clump. Harries] (120001) studied 
the time-variability of a clumpy wind, again with a Monte-Carlo 
code, but did not investigate the effect of changing stellar/wind 
properties, nor the polarimetric variability on longer timescales. 
Both of these studies concluded that the clumps must be opti- 
cally thin, as increasing the clump optical depth begins to have 
a conspicuo us effect on tim e-resolved observations of spectral 
line-profiles. iLi etal (2000) used an analytical model to investi- 
gate the polarimetric variability of WR s, but again did not study 
the effect of different stellar parameters. Richards on et al.l (1996) 
studied specifically the effect of a clumpy wind on the ratio 
of photomet r ic to p olarimetric variability in the case of WRs. 
Bro wn et"aT] (Hp00) investigated the 'redistribution' of a smooth 
outflow into clumps, simulating the pile-up of material caused 
by e.g. shocks in the wind. They concluded that small-scale re- 
distribution close to the star could not produce significant lev- 
els of polarization. This implies that clumps producing short- 
timescale polarimetric variability, which are necessarily close to 
the star, must be localized over-densities of ejecta, as opposed 
to the result of a smooth outflow being 'snow-ploughed' by a 
shock-front. 

In this paper we use an analytic model to build on these pre- 
vious studies and fully investigate the polarimetric variability 
produced by a clumpy wind across a range of stellar parameters 
typical of massive evolved stars. We place particular emphasis 
on LBVs - the group of stars which display the strongest polari- 
metric variability. A full description of the model is presented 
in Sect. [2] along with a few basic applications to geometries ex- 
plored in other work to serve as validation. The relevant input 
parameters are discussed in Sect. [3] leading to the model results 
and discussion in Sects. |4]and[5] 



2. 1 . Model geometry 

The c oordinate system used is that defined in lBrown & McLean] 
d 1977b . and is illustrated in Fig.Q] The angle of the plane of scat- 
tering to the observer ;f can be described in terms of the polar (#) 
and azimuthal (<p) angles of the star's reference frame, and the 
inclination of this frame to the plane of the sky i: 

co&x - - cos cos / + sin 6 sin i cos (1) 

The position angle (PA) of the polarization, , is described by, 

cos 9 sin i + sin 6 cos i cos d> 



sin 



{(cos 9 sin i + sin 6 cos / cos <f>) 2 + sin 9 sin <p] 



0.5 



(2) 




-X 1 



Fig. 1. The coordinate sy stem used in this model, as defined in 
Brown & McL eanl dl977b . A beam with direction (9, (p) scatters 
through an angle x towards the observer. The observer's refer- 
ence frame is rotated through an angle / with respect to the star's 
reference frame. 



2.2. Polarization of axi-symmetric material 

From iBrown et al .1 (2000), which uses the expr essions of 
IBrown & McLeanl d 19771) and lCassinelli et al.l dl987b . the polar- 



ization produced by electron scattering material in a stellar wind 
as a function of distance to the star r, polar angle 9, and az- 
imuthal angle (p, is described by, 



2. Description of the model 

In this analytic study we f o llow the approach o f, amongs t others , 
IBrown & McLeanl (U977I) : IBrown et al.l (120001) ; ILi et al.l §000). 
The strategy for the construction of the model is as follows: the 
polarization due to a single clump is determined as a function of 
distance from the illuminating star. Through the velocity-law of 
the wind, an expression for the polarization as a function of time 
is found. Multiple clumps are added, and the total polarization 
found from the vector sum of the polarization of each individual 
clump. The clumps are then given a typical ejection timescale, 
and the polarization as a function of time is calculated as new 
clumps are ejected and move through the wind. After the simu- 
lation reaches a steady-state, the polarization is monitored for a 
further period of time to determine the typical polarization level 
and its variability. 



3 r^ 2 

P(fi, r) = — o" T sin 2 ^ I (1 - 3j/)dfi 



x J n e (r)-yi -(-^-j dr 

where cr T is the Thompson scattering cross-section, x is the an- 
gle between the plane of scattering and the observer, fi = cos 9, 
R+ is the radius of the star, and n e is the electron number density 
which in this case is a function of r only. The square-root factor 
within the r integral is known as the finite-star correction-factor 
(Cassinell i et al.ll 19871) . and accounts for the finite stellar 'disk' 
seen by the scattering material. This formalism can be used to 
determine the polarization produced by different scattering ge- 
ometries, and an illustration of this can be seen in Fig. [2] 
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Fig. 2. The geometry assumed in Eq. (|3}. The electron density 
is assumed to be symmetric about the central axis (i.e. over all 
(/>), such that by specifying the limits r\ to r 2 and 9\ to 62 we are 
integrating over the volume swept out by the shaded area, in this 
case a 'shell'. 

2.3. Clumps 

The amount of polarization a clump produces is determined by 
how much light is intercepted, as well as the clump's size and 
morphology. Arbitrary changes in clump morphology can pro- 
duce minor changes in P if the clump is close to the star, due to 
the effect of the finite size of the illuminating star. This effect is 
small, however, and it is a reasonable approximation to say that, 
for a constant column-density, halving a clump's size (and hence 
doub ling its density) pr oduces no difference in polarization (see 
also iBrown et alJfeOOOl) . For this reason we concentrate on the 
simplest morphology for the clumps, that of constant angular 
size and thickness. In a radiatively driven wind this scenario is 
plausible, as each particle in the wind is driven radially outwards, 
and any inhomogeneity will grow in size with the square of the 
distance. 



2.3.1 . Polarization of a single 'shell' clump 

The geometry of Fig. [2] and Eq. (0 are used as a starting point. 
The angular size remains constant with distance from the star, 
so \x can be integrated from fii = to fi2 = 1 - AQ./2n, where 
AO. is the angular size of the clump. The electron number den- 
sity is parameterized in order to conserve mass for a clump with 
constant angular size, thickness and number of electrons, i.e. 



(3) 



where «o m the initial number density of the clump at r = R+. 

Substituting Eq. © into Eq. O) and setting x = r/R*, the 
polarization of a single clump at a distance x' from a star is, 



Id 



where Ar is the dimensionless thickness of the shell. 

As this is an analytical study, multiple scattering is not taken 
into account in this equation. When multiple scattering effects 
start to become relevant, i.e. when r e > 1, the results from 
Eq. would become unphysical. To assess the importance of 
multiple scattering, we need to monitor the optical depth, and we 
re-write Eq. © in terms of the initial optical depth per clump tq, 
i.e. when r = R+: 



TO 



3 n 

16 T J n 



(r)dr 



—o- T n Q R+K 
10 



where, 



dx 



1 



1 + Ar 



=s Ar (if Ar 1) 

In order to relate the initial optical depth to other observ- 
ables, the initial optical depth per clump can be expressed in 
terms of the stellar parameters. The mass-loss rate and clump 
ejection rate, N, are related via, 



M = NN e /j e m H 



(4) 



where N e is the number of electrons in each clump, yu e is the 
mean mass-per-particle and m# is the mass of a proton. To pa- 
rameterize the ejection rate in the general case, the formalism of 
iLi et al.l (l2000h is adopted by defining the ejection rate per wind 
flow-time N, where ta = R^/v^: 



N = Ntfi = NR+fVoa 

Therefore, the number of electrons in each clump: 

MR* 
N e - * 

fi e m H Nvoo 
The initial volume per clump is, 

V d = AORlAr 



(5) 



(6) 



(7) 



where AO is the solid angle 2x(l - fj.2), and Ar is the clump 
thickness in units of and is assumed to be <scl. The electron 
density per clump is then, 



N, 



M 



Equation can now be re-written as, 



n T ■ 2 / 3x 

P = T sln ^(y"2-^ 2 ) 
Ar 



(8) 
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Fig. 3. The polarization as a function of distance for a single 
clump for three values of t. The clump size is AO. = 0.04, Ar = 
0.01; and sin 2 ^ = 1. The polarization falls-off as roughly l/r, 
with the attenuation at small r due to the finite size of the illu- 
minating star. It can be seen that at distances of r/R+ ~ 5 the 
polarization has already decreased from the maximum value by 
a factor of ~10. 



where, 

TQ_ _ 

Ar ~ ' 



M 



0~T 

16 p e m H A/v c 



ACIR+Ar 



(9) 



The results of the polarization of a single clump as a function of 
distance are shown for different optical depth values r in Fig. [3] 



2.3.2. Polarization of a single clump as a function of time 

Having determined the polarization as a function of distance 
from the central star, we calculate the time behaviour of the 
polarization. The velocity of the wind as a function of distance 
from the star is described by the so-called 'beta-law', 



v(r) = Voo 1 - 



bR-k 



(10) 



where v m is the terminal velocity of the wind, b is a constant 
which dictates the initial velocity of the wind vo (i.e. at r = 7?*), 
and /3 is a constant which determines the wind's acceleration. 
Hot supergiants are typically thought to have 1 < B < 2 (see e.g. 
iHerrero et a l. 2002), so for ease of calculation a value of ft = 1 
is used throughout most of this study. By rearranging Eq. ( ITOl ). 
setting p = 1, and substituting x = r/R+, we find an expression 
for the time since ejection t as a function of distance x and the 
velocity-law parameters Vco,/3, and b, 



J v(r) Vco 



[x + b ln(x - b) + c] 



(11) 



where c is a constant of integration, which is determined by forc- 
ing x = 1 at t = 0. The polarization as a function of time can now 
be calculated by solving Eqs. and (fTTT i parametrically. The 
initial velocity Vo is set to the isothermal sound speed at the tem- 
peratures considered. Altering vo by a factor of two was found 
to have little difference on the results of a given simulation, so 
for simplicity the value vq = 10 km s _1 is used throughout this 
study. 




10 



20 



30 



Fig. 4. The polarization as a function of time (in units of wind 
flow-time) for three different values of p. Top: for slow termi- 
nal wind speed; Bottom: fast terminal wind speed. Same clump 
parameters are used as for Fig. [3] 



The effect of increasing yS is shown in Fig. [4] The plot shows 
the polarization as a function of time for three values of and 
two separate ratios of vo/v'co. As increasing B increases the time 
the clump spends in the inner regions of the wind (the 'dwell 
time'), the polarization remains high for longer. The effect is 
more pronounced when vo/voo is higher. 

2.3.3. Polarization of multiple clumps 

Clumps are ejected radially from a random location on the stellar 
surface, with < fyj < 2n and - 1 < cos 9j < 1 for the /'th clump. 
For a given inclination angle, the polarization and its position 
angle due to a single clump as a function of time are calculated 
through Eqs. ([T), (|2), (O, and ( fTTT i. Clumps are ejected over a 
typical timescale T e j, such that T e j = R^/Nv^. The effect of 
randomizing 7Vj was experimented with, such that Tj = T ej + 
cr e j. This was found to have no effect on the results, and was 
computationally more expensive as it required finer grids to be 
set up. Consequently, a uniform T e j was used throughout each 
simulation. 

To calculate the total polarization due to all clumps at any 
one time, the polarization due to each clump must be converted 
into the linear Stokes parameters, 



Q = Pcos2© 
U = Psin2© 



(12a) 
(12b) 
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where © is perpendicular to the plane of scattering. The total 
polarization at any one time is then determined from the sum of 
the Stokes vectors of each of the N clumps in the wind, 



N 



©tot = - tan 1 



I N 



V 7=1 



I N 



./=! 



(13a) 



(13b) 



In order to determine the typical polarization of a model run, 
the simulation is run until it reaches an equilibrium, and the 
first clumps to be ejected no longer contribute significantly to 
the overall polarization. This occurs when the first clumps reach 
r/R-,, ^ 6 (see Fig. [3). The polarization is then still monitored 
for a further 30 flow-times, and the time-averaged polarization 
(P) over this period, together with the r.m.s. variability AP, are 
recorded. As the simulation is stochastic by design, it can some- 
times throw up atypical results. To compensate for this, the pro- 
cess is repeated several times for each set of parameters, and 
the mean is taken. The experimental uncertainties in (P) and AP 
were taken to be cr/ V", where <x is the r.m.s. scatter of the results 
of n simulations. 



2.3.4. High ejection-rates 

The goal of this work is to reproduce the polarimetric variability 
of hot stars - both on long timescales, as found by the repeated 
meas urements in th e lite rature; and shor t -times cales, as found by 
e.g. lHavesI dl984al) and ISt.-Louis et alj (119871) . At low ejection 
rates, individual ejection events may be resolvable by nightly 
monitoring, and the model would have to be dynamic to simulate 
this, i.e. follow the polarimetric behaviour in real-time as clumps 
move through the wind, whilst new clumps are ejected. 

However, at high ejection rates the distribution of clumps 
close to the photosphere will look completely different from 
night-to-night, and so will the polarization. In this situation, 
there is nothing to be gained from computationally-expensive 
dynamical modelling as ejection events are unresolvable. In the 
high-A/ regime, calculations are therefore performed using a 
static model. Here, as the polarization at each observation is un- 
related to the last, a series of wind 'snap-shots' are generated, 
which represent a sample of independent observations. Clump 
distances are randomly-generated according to the ejection rate 
and velocity law, and the total polarization is calculated. This is 
repeated 10 3 times, and P and AP calculated. 

Whilst being unable to predict the time-evolution of individ- 
ual clumps, this static method has the advantage of being able 
to simulate a much broader range of N values - at drastically 
reduced computation times. 



2.4. Optically-thick clumps 

Equation ([3]) shows that, in this study, the amount of polarization 
produced by a clump is proportional to the clump's optical depth, 
as denser clumps scatter more light. 

However, this will only hold as long as the clump is optically 
thin. When the optical depth becomes larger than ~ 1, multi- 
ple scattering becomes important. In this regime, the relation in 
Eq. (01 breaks down as a photon which is scattered twice cannot 
produce twice as much polarization. 

To investigate the effect of optically-thick clumps, a first- 
order approximation is made from the Monte-Carlo work of 



iRodrigues & M agalhaes (2000). They studied the polarization 
produced by starlight scattering off clumps with optical depths 
as high as r e = 10. Their results show that, in the low-T e regime, 
the polarization rises almost linearly with optical depth, until t 
reaches ~ 1 . At high-r e the polarization flattens off. This can be 
understood in terms of multiply-scattered photons 'forgetting' 
where they came from, with only the singly-scattered photons 
carrying the scattering geometry information. This explains the 
slight fall-off of P with high r, as less and less of the observed 
photons are singly-scattered. 

We approximate the results of Rodrigues & Magalhaes 
(2000) such that polarization is linear with optical depth up to 
T = 1, and then constant for higher r. In this situation, as an 
optically-thick clump moves away from the star its density de- 
creases as 1 /r 2 , but the polarization remains constant while the 
optical depth is larger than 1. When the clump reaches a dis- 
tance such that t - 1, its polarization decreases rapidly in the 
same manner as an optically-thin clump. 

In making t his a pprox imation, we note tha t both 
ICode & Whitoevl d 19951) and IRodrigues & Magalhaes! d2000h 
considered the case of spherical clumps, while in this work 
we consider only flattened, 'shell' clumps. Although we argue 
in Sect. 12.31 that clump shape is a redundant parameter in the 
optically-thin limit, the relative invariance of clump morphol- 
ogy to polarization may not hold for optically-thick clumps. 
For example, in the case of a spherical clump, the limb will 
always be optically-thin to some extent, whereas the optical 
depth of a 'shell' clump will always be the same along all 
radial paths from the star. This problem bears some similarity 
to the challen ges faced when modelling clouds in planetary 
atmospheres (Ivan de Hulsllll980l) . The impact and limitations of 
our approximation in the context of the results will be discussed 
further in Sect. 1441 



3. Model parameters 

For a given mass-loss rate and ejection timescale, the key input 
parameters are: the terminal wind velocity, the clump size and 
optical depth, the stellar radius, and the ejection timescale. Each 
of these parameters are discussed below. 



3.1. Clump size 

The size of each clump is determined by its angular size AO. and 
its thickness Ar. Provided Ar is small, it has no effect on the 
polarization per clump as the optical depth remains the same. 
Varying AQ. has a small effect, as for large A£2 the polarization 
at the edge of a shell slightly cancels that of the opposite side. 
This can be seen in the simulations (Sect. l4~2~2l . 



3.2. Terminal wind velocity Vco, and stellar radius P* 

The polarization produced by a single clump falls off quickly 
with distance to the star (see Fig. [3). The amount of time that 
each clump will spend close to the star (the 'dwell time') is there- 
fore a major factor in the overall level of polarization. The dwell 
time is affected not only by speed of the outflow, but also the size 
of the star. This is therefore reflected in the star's wind flow-time 

ffl = P*/Voo- 

The terminal velocity is typically a few times the stellar es- 
cape velocity, which in turn depends on the effective gravity of 
the star: 



eif - gN - gmd 



(14) 
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where is the Newtonian gravity GM/R 2 and g la d is the out- 
ward radiative acceleration due to the continuum opacity of free 
electrons in the ionised stellar atmosphere, 



(15) 



4-7TK+C 

The effective mass is then defined as, 

M eff = (l-F e )M* (16) 
where F e = grad/gN- The terminal wind velocity is now, 



fiBSjV'e 



= Q 



BS.I 



2GM, 



eff 



R* 



(17) 



where Qbsj is the "bi-stability jump" ratio whi ch is pre- 
dicted to depend on the star's effective te mperature dVink et al.l 
19991). According to empirical studies (lLamers et al.l 11995; 
Kudrit zki & Pulsll2000h . <2bsj takes two distinct values, 



gBSj-1-3 (r eff ^ 21000K) 
£bsj - 2.5 (r eff £ 21000K) 

As the stellar radius and effective temperature are related via, 

U = 4nRlcr B T^ (18) 

the terminal wind velocity can be found from the star's mass, lu- 
minosity and temperature. In the case of LBVs, this allows us to 
predict how the terminal wind-speed varies over the S Doradus 
phases which are thought to occur at constant bolometric lumi- 
nosity. 

We note that a recent study of blue supergiants over the bi- 
stability regime indicates that 2bsj may not be a step-function, 
but in stead it may decrease smoothly with decreasing temper- 
ature (Crowth er et al.l 120061) . The binning of the data by effec- 
tive temperature by Lamers et al. may give the impression that 
2bsj behaves like a step function. One of the main reasons Qbsj 
may not be a step function is that one is comparing objects with 
very different stellar parameters. For individual objects, as is 
relevant for our study, the wind velocity may nonetheless be 
expected to 'flip-flop', especially in the doma in of LBVs (e.g. 
iPauldrach & Pulslll990tlVirik & de Koterll2002l) . 

4. Results 

In the first part of this sectio n, we compare our model to the re- 
sults of a similar study by iLi et ail (120001) . Next, the effects of 
the various input parameters are discussed one-by-one. By con- 
straining various stellar parameters to those of the well-studied, 
polarimetrically-variable P Cyg, the more detailed parameter- 
space is explored. 

4.1. Comparison with previous work 

ILi et al.l {2000) investigated the polarimetric variability of Wolf- 
Rayet stars (WRs). Figure [5] shows a comparison between their 
results and the results of our model. Qualitatively the two stud- 
ies show a similar trend of steadily-falling (P) with increasing 
Af, behaviour also found by Richardso n et al.l (U996). However, 



quantitatively the two differ by a factor of ~5. The reason for this 
is their different treatment of the clump density. 

Li et al. fixed the clump angular size AO., clump thickness 
Ar, and the electrons per clump A^ e at constant values. Therefore, 
the volume of the clump V c \ should increase with the square of 
the distance, 



Vol oc r 1 



(19) 



They also stated that, following mass continuity, they scaled 
their electron density as 



H e a l/(r v) 



(20) 



These statements are not consistent: if the number of electrons 
per clump N e = n s V c \, then from Eq. d20t . 



N e oc l/(«eV) 



(21) 



Clearly, the number of electrons cannot be both constant and 
proportional to 1/v. For the system of equations to be self- 
consistent, the clumps must either, 

i) be allowed to grow in thickness as they propagate through 
the wind, due to the differential acceleration between the 
front and back of the clumps; or 

ii) have a fixed thickness and an internal clump density which 
scales as 1/r 2 . The space density of the clumps still falls as 
l/r 2 v 

Li et al. integrated from (x\ = x') — » (x2 = x'+Ar), where Ar 
was constant, whilst allowing the electron density to fall as 1/r 2 v. 
They therefore only integrated over the front Ar of the clump, 
systematically underestimating the polarization per clump. This 
explains their lower levels of (P) compared to ours. 

The differences in polarization per clump as a function of 
distance from the star is illustrated in Fig. [6] Three scenarios are 
plotted: 

1 Clump density oc 1/r 2 , with fixed AO., Ar. 

2 Clump density oc 1 / r 2 v, with fixed AQ, Ar. 

3 Clump density oc l/r 2 v, with fixed AQ, but with Ar growing 
proportional to v. 

It can be seen that there is a significant difference be- 
tween scenario (2) and the self-consistent scenarios (1) and (3). 
Quantitatively, the difference between the results of Li et al. and 
this work can be understood in terms of the maximum polar- 
ization per clump P max . In the limit where the ejection rate is 
much shorter than a flow-time, i.e. N » 1, the time-averaged 
polarization (P) is proportional to P max . The maximum polar- 
ization per clump in scenario (2) is ~0.1%, while for (1) and (3) 
Pmax ~ 0.5%. This difference could explain the factor of ~ 5 
difference between this work and Li et al. for this set of model 
parameters. 

It can also be seen from Fig. [6] that polarimetrically there 
is little difference between scenarios (1) and (3). This is to be 
expected, as both scenarios integrate over the same amount of 
material and have the same optical depth. The minor difference 
when x is small is due to the finite star correction term. The 
physical differences between scenarios (1) and (3) are as fol- 
lows: (1) represents a thin shell of material being lifted off the 
surface of the star; while (3) describes a localised density over- 
enhancement, the material of which moves through the wind in 
the same manner as the ambient material. We note that the two 
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Fig. 5. Comparison between the results of iLi et alj (12000b and 
the results of this work when divided by a scaling factor of 5. 
The top panel shows the time-averaged polarization as a function 
of clump ejections per wind flow-time, whilst the bottom panel 
shows the variability in P. The results of the present model have 
been scaled to highlight the fact that while the two models agree 
qualitatively, quantitatively they differ. Parameters used are Voc = 
1800km s _1 , R* = 10R@, b = 0.995, and \og(MIM & yr l ) = -4. 



scenarios have very different implications for the clump filling- 
factor f c \\ (3) has / C ] constant at all radii; while (1) has f c \ which 
grows with distance from the star to some asymptotic value. 
Figure [6] shows that while the two scenarios differ in the phys- 
ical situation they describe, the resulting polarization is almost 
identical. As the computation times for (1) are much shorter, and 
integrals easier, this scenario is preferred throughout this study. 

4.2. Effect of various input parameters 

The free parameters in the model are the clump angular size AQ., 
ejection rate N, and the stellar parameters fi e , R+, M, v M , and Vo 
(i.e. b). Of these, AQ, fi e and M will be fixed throughout the pa- 
rameter study, so the effect of these variables is discussed below. 

4.2.1. Mass-loss rate 

The effect of varying M for a fixed set of model parameters 
is shown in Fig. [7] In the optically-thin limit, the polarization 
per clump is directly proportional to M (Eqs.|9]and|9]l. Figure|7] 
shows, as expected, a linear relationship between M and (P) in 
the optically-thin regime, and a linear least-squares-fit is over- 
plotte d to il l ustrat e this. This linear relationship was also found 
bv lLi et all d2000b . As the mass-loss rate increases the clumps 
become optically-thick, therefore we would expect the trend to 
flatten-off for extremely-high mass-loss rates. 

Also shown in Fig.|7]is the calculated polarimetric variabil- 
ity, AP. This is also linear with M, with AP/(P) consistently 
found to be ~ 0.5 — > 0.6, not just as a function of mass-loss rate 
but for all simulations. This can be understood as follows: if P 
is truly random, one would expect it to have the appearance of 
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Fig. 6. Comparison between the three scenarios described in the 
text for calculating polarization as a function of distance for a 
single clump. Scenarios (1) and (3) are self-consistent and pro- 
duce almost identical results, while (2) is different. Curve (2) 
peaks at ~0.1% extremely close to the star, while (1) and (3) 
peak at ~0.5% at a distance of 1.2/?*. It is speculated that it is 
this factor of 5 difference in maximum polarization per clump 
which is responsible for the differing model results (see Fig.[5j. 
The curves are evaluated for t — 1, sin 2 ^f = 1, AQ = 0.04, and 
b = 0.995. 
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Fig. 7. The effect of varying mass-loss rate on the time-averaged 
polarization and its variability in the optically-thin regime. The 
model parameters used are Voo = 180km s~', R+ = IOORq, AQ 
0.04. Ar 0.01 



noise, with the amplitude — > P m!ix to be roughly 3cr. The time- 
averaged value of P will then be between these two at ~ 1 .5cr. 
Therefore, the quantity AP/(P) ~ 1/1.5 ~ 0.67. 



4.2.2. Clump size 

If the number of scatterers in a clump is conserved, then in 
theory changing the clump's size should have little effect, as 
the same amount of light will be scattered. However, there is 
a slight trend of decreasing polarization for larger clumps: con- 
sider a 'shell' clump ejected perpendicular to the observer's line- 
of-sight. When viewed on the plane of the sky, light scattering 
towards the observer from the two opposing limbs of the clump 
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Fig. 8. The effect of varying clump angular size on the time- 
averaged polarization. Overplotted is the relation in Eq. ( 1221 ). to 
which (P) is expected to be proportional (dashed line). Errors on 
the data points are from the uncertainties in the means of seven 
independent simulations. Input parameters are the same as used 
in Figs.[5]&|7] 



will produce polarization with slightly different position-angles. 
For a clump with a large angular diameter, the position-angles 
of these two beams will tend toward being perpendicular to each 
other, thus we expect less net polarization from larger clumps. 
Quantitatively, for otherwise constant input parameters, from in- 
spection of Eqs. (0 and (0 one would expect: 



(P) cc (j^-^/Afi 



(22) 



where fi2 = 1 - AQ/27T. The effect of varying AO is shown in 
Fig. [8] with the relation in Eq. d22l scaled and overplotted. The 
behaviour of (P) follows this relation well, with the scatter due 
to the stochastic nature of the process. 

Large values of AQ (^0.5) are likely to be physically unreal- 
istic, and the fractional change in (P) between AQ = 0.01 — > 0.5 
is small (Fig. [8}. There fore, we fix the angular size of the clumps 
to the value chosen bv lLietal.1 d2000l) . that of AQ = 0.04. 

4.2.3. Mean mass-per electron 

The process considered in this study is electron-scattering, how- 
ever the electrons in each clump contribute a negligible amount 
to the clump's mass. It is assumed that there are nuclei to account 
for each of these free electrons, e.g. if the wind material was 
pure H then the mass of a clump would be N^mu- However, the 
material in the winds of massive evolved stars will have a non- 
negligible content of heavier elements, so the mass per clump is 
determined by the mean mass-per-electron, /i e . From Eq. (0 it 
can be seen that the optical depth (and hence the polarization per 
clump) is inversely proportional to fi e . The reason for this is that 
greater values of yu e mean fewer electrons in a clump, reducing 
the electron-scattering optical depth. 

The elemental abundances in a star's wind are determined 
to first order by the chemistry of the natal ISM and the evo- 
lutionary phase of the star. From recent abundance studies of 
hot supergiants it is assumed that the elements heavier than he- 
lium contribute a n egligible amount to a clump's mass (e.g. 
ICrowther et alj 120061) . Determinations of the helium number 
fraction Y - «He/«H ar e not straight-forward, due to subtle 
dependences of He line-strengths on the details of radiative 



transfer calculations (see e.g. iMcErlean et al.lll998l). Studies 
of LBVs have yielded v alues of Y - 0.3 - 0.5, (ICrowfherl 
ll997tlNaiarro et aLlll997l) ; while more recent studies of OB su- 
pergiants with upg raded models suggest slightly lower values 
(Y = 0.1 - 0.3, iHerrero et al.ll2002t iRepolust et al.ll2005l) . This 
implies a value of yu e ~ 1.3 - 1.6. As the uncertainty in the mean 
mass-per-electron is small compared to that of other parameters, 
a constant value of yu e = 1 .5 is adopted throughout the study. 

4.3. Polarization as a function of stellar radius 

Following the above arguments, this section investigates the po- 
larization produced by a clumpy wind across the range of stel- 
lar radii inferred from LBV temperatures. For a given mass and 
constant bolometric luminosity, the stellar radius and therefore 
effective temperature and terminal velocity can be varied. Then, 
for a given mass-loss rate, the clump ejection timescale and opti- 
cal depth can be varied, and the resulting polarization calculated. 
Values of L = 10 6 L Q , M = 3OM , and log(M /Af yr-') = -4.5 
are us ed, which are typical of the LBV AG Car dVink & de Koterl 
2002). The effect of varying M is significant, as it is proportional 
to the time-averaged polarization (see previous section). L and 
M, in this parameterization, only affect the terminal wind veloc- 
ity, and so are redundant parameters if Vco is known. The polar- 
ization as a function of v M , and T e ff for a range of ejection 
timescales N is shown in Figure|9] 

4.3.1. Varying ejection timescale 

If the mass-loss rate is conserved, increasing the number of 
clumps ejected per wind flow-time N will decrease the optical 
depth per clump. For low N, the wind will consist of a small 
number of dense clumps. For increasing N, the wind consists 
of a larger number of clumps with smaller optical depths, until 
eventually the wind tends toward a smooth, spherically symmet- 
ric outflow. Figure [9] illustrates this. It can be seen that while 
the level of polarization decreases with increasing N, even at a 
high ejection rate there is still a detectable level of polarization 
(>0. 1 %) for the wind parameters used. 

4.3.2. Varying clump optical depth 

For constant clump size, decreasing the initial optical depth per 
clump will also decrease the mass per clump. Therefore if the 
mass-loss rate is to be conserved the ejection rate must increase. 
The wind will be highly clumped at To = 1, and will tend towards 
a smooth outflow at low values of r. For the model parameters 
of Fig.0 N is linearly related to t via, 



■ = 3.6X10** (r^rfeV 
Ikms-'/ \R j 



(23) 



Evaluating this expression for the lowest curves in Fig. [9] optical 
depths are in the region tq = 0.05 — > 0.1. As Voo oc 1/ y7?7 (Eq. 
[TTI i, optical depths are highest for small radius, high terminal 
wind-speed stars. 

4.3.3. Overall behaviour of P with tf* 

It can be seen from Fig.0that, for a given ejection-rate per wind 
flow-time there is a general trend of increasing polarization with 
decreasing radius. This can be understood as follows: as the size 
of the star is decreased, the clumps become smaller and denser, 
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Fig. 9. Time averaged polarization for constant values of TV, as 
a function of terminal wind speed (top), /?* (middle), and effec- 
tive temperature (bottom). The stellar parameters are the same as 
those used throughout the parameter study, namely L = 10 6 L Q , 
M = 3OM , and logCM/Moyr 1 ) = -4.5. 



producing more polarization per clump. Around the bi-stability 
zone, the terminal velocity increases, and clumps spend less time 
close to the star: the overall polarization drops again. 

The overall levels of polarization produced in the simulations 
are high - even at high TV ~ 300 where the wind is tending to- 
ward a smooth outflow, observable residual levels of polarization 
(~0.3%) are produced. This is discussed further in Sect. [5] 



4.4. Optically-thick clumps 

In this and the following section, the parameter-space is explored 
in greater depth. To do this, the mass-loss rate, stellar radius and 
terminal wind-speed are fixed at those of the LBV P Cyg, the sta r 
for which most pol arime tric observations exist. lHavesI (119851) . 
iTavlor et ail d!991l) and iNordsieck et alj d2001l) each present 
many individual measurements, finding (P) a 0.3 ± 0.15%. 
Indeed, it is an extremely well studied object generally, with 
many independent measurements of stellar / wind parameters in 
the literature. As the star hasn't been significantly variable for 
^200 years, this means that the input parameters of stellar ra- 
dius, mass-loss rate and terminal wind speed can be constrained. 
This then allows further exploration of other details which may 
or may not be significant. 
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Fig. 10. Time-averaged polarization over a broad range of TV. At 
N ~ 20, the optical depth per clump exceeds unity, and the over- 
all polarization fall s (see text for details ). The P Cyg wind pa- 
rameters derived by Naia rro et al. (1997) are used. The observed 
level of polarization for P Cyg is marked by the dash-dotted line. 



Th e parameters used are those derived by iNaiarro et alj 
(119971) . namely 

R*/R Q = 75, Vco = 18510ns- 1 , logCM/Moyr 1 ) = -4.52 

These values were found through spectroscopic analysis with 
non-LTE modelling. The model took no account of wind- 
clumping, therefore the mass-loss rate may be an overestimate. 
Also, the authors note that there is some degeneracy between 

v M and M such that M/(VooR+ 5 ) ~ constant, but constrain 
R+ through other derived values in the literature. 

Figure [10] shows the time-averaged polarization over a 
broader range of TV than in Fig|9] At high TV, the static model 
was used to calculate P, and it can be seen that the level of po- 
larization observed for P Cyg is reached at TV ~ 2000. 

At low TV, the density -per-clump becomes very high in order 
to conserve the mass-loss rate, and exceeds unity below TV ~ 20. 
At this point, the model switches to the optic ally- thick approx- 
imation, that is the polarization per clump plateaus at r > 1. 
Therefore at TV < 20 the overall polarization begins to fall again. 
We again point out that, in the optically-thick regime, the polar- 
ization may become sensitive to clump morphology, and hence 
the optically-thick part of the curve in Fig. [10] (left of TV ^ 10) 
may not be as quantitatively robust as the thin-clump regime. 
However, the qualitative behaviour of a turnover in (P) is in- 
evitable, as there is only a maximum amount of polarization a 
single clump can produce. Therefore, we expect that as clump 
density increases and the number of clumps becomes smaller, 
there will be two solutions to match the observed level of polar- 
ization - the optically-thin, and optically-thick cases. 

Extrapolating from Fig. [TO] the optically-thick solution for 
P Cyg is reached at at TV S5 0.01. As the flow-time for P Cyg is 
about 3 days, this would imply ejection rates of only a few per 
year even allowing for a factor of 10 difference in the calcula- 
tions. This is deemed to be unlikely, for reasons which will be 
discussed in greater detail in Sect. [5] and we therefore favour the 
optically-thin solution. 
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Fig. 11. The polarization produced when the wind acceleration 
parameter/? is increased, lowering the acceleration of the wind. 
Stellar parameters are the same as Fig.fTOl 



4.5. Lower wind acceleration 

In Fig. [TT] the effect of increasing the wind acceleration param- 
eter fj is explored. It can be seen that, in the case of P Cyg, 
increasing f3 makes very little difference. This is related to the 
ratio vo/voo: from Fig. [3] it is clear that the greatest amount of 
polarization is produced when the clump is within ~ 6R* of the 
star. Figure [4] shows that when vo/v^ is relatively large, altering 
fj makes virtually no difference to the dwell-time. Only when 
vo/vco is small is the dwell-time significantly affected by fj. It 
is noted that at larger v^, increasing fj has a more pronounced 
effect, as would be expect from the low vo/vco case shown in 
Fig-El At vo/vco ~ 1/180, the polarization is increased by a fac- 
tor of ~ 1.5 when fj is increased to 2. 



5. Discussion 

5.1. A few massive clumps or many tiny clumps? 

In the approximation of a wind consisting of uniform clumps 
with no inter-clump medium, the densities of the individual 
clumps need not be large to produce the levels of polarization 
observed in massive stars. Polarimetric variability on the scale 
observed in LBVs can in principle be reproduced by two distinct 
regimes of N: either ^0.1 in the optically-thick case, or ^1000 
in the optically-thin case (see Fig.fTOt. 

In the former, the time between ejections would be of the 
order of months (the flow-time of P Cyg is about 3 days ), and 
one w ould expect ejection events to be seen in the data of lHavesi 
( 1985) whose temporal resolution was of order 1 day. Also, the 
density of the clumps would be so large that the recombina- 
tion emission would almost certainly be conspicuous in time- 
resolved spectroscopy (quantitative analysis of this is beyond the 
scope of the present work). For this reason, the higher value of 
TV is considered to be the most likely value of the two. 

At the inferred ejection rates of N ^ 10 3 , the inner wind con- 
sists of thousands of low density clumps. In this regime, one may 
expect the clumps to cancel each other out and leave a net po- 
larization of zero. However, only a slight imbalance is required 
to produce residual polarization. For example, at TV = 2000 and 
the stellar parameters of P Cyg, the maximum polarization per 



clump is 0.006%, and there are ~ 7 x 10 3 clumps in the in- 
ner 2R*. A Poissonian l<x deviation from spherical symmetry 
of ~ 80 clumps is enough to produce a residual polarization of 
{P)= 0.4%, which is the observed level of P Cyg. Therefore, 
the polarization may simply result from the statistical deviations 
from spherical symmetry: a slightly -fragmented wind model. 

If observable levels of polarization are produced when a 
wind is so close to being homogeneous, it is reasonable to ask 
wh y intrinsic pol arization is not present in all the stars observed 
by dDavies et alj|2005l) . The answer may be related to the star's 
mass-loss rate: as shown in Fig. [7] the polarization produced by a 
given simulation is proportional to M if all other parameters are 
fixed. Therefore, if a star's ejection rate were fixed but its mass- 
loss rate were reduced from 10~ 45 to 10~ 6 Moyr _1 , its average 
polarization could go from being as large as any observed to be- 
ing completely undetectable. This suggests that intrinsic polar- 
ization may be ubiquitous among hot stars with mass-loss rates 
greater than 10 M yr _1 , as a measurable level of polarization 
can be produced from a statistical deviation from spherical sym- 
metry. 

Th is conclusion is supported by the findings of D avies et~ai] 
(2005). Intrinsic polarization is more likely to be found in stars 
with strong Ha emission; and for LBVs, whose terminal wind 
velocities are roughly the same, Ha equivalent width may be 
considered to be a crude measure of mass-loss rate. It is also con- 
sistent with the lower inci dence of polarization in O- supergiants 
(Osgs) and WRs, found in lHarries etaT] dl998ll2002l) Osgs have 
typical mass-loss rates of 10 _5 M G yr _1 dPuls et aLl u006): and 
while WRs are considered to have M ~ 10~ 5 — > 10 _4 M Q yr _1 , 
ironically this may be redu ced due to the effect of clumping 
dHamann & KoesterkdH99 8). The impact of clumping on LBV 
mass-loss rates is yet to be studied with the latest versions of 
model atmospheres. 

Of course, there is the possibility that the wind is a combina- 
tion of the two scenarios, where a slightly fragmented wind pro- 
duces low-level variability, combined with sporadic ejections of 
optically-thick clumps. Such a wind could be simulated by re- 
placing the uniform clumps of the present model with clumps 
whose optical depths are determined by a power-law distribu- 
tion. This would involve the introduction of an extra free pa- 
rameter, the power-law index, and so it would be reasonable 
to look for observational evidence for such a scenario first be- 
fore increasing the degenerate parameter space. The polarimet- 
ric signature of a dominant dense clump would be a jump in 
polarization and a PA which remained constant for several flow- 
times (few weeks for P Cyg), until its effect is 'washed-out' 
by the smaller clump s close to the star. Also, it was shown by 
Igna ce & Churchweli ([2004) that hierarchical clumping of ion- 
ized outflows can lead to atypical spectral indices in IR and 
radio free-free emission. Such spectra have slopes intermedi- 
ate between the optic ally- thick/thin limits over several orders of 
magnitude in frequency. Hence, a combination of polarimetric 
monitoring and dynamic modelling, plus continuum radio obser- 
vations, may be used to investigate hierarchical wind-clumping 
further. 



5.2. Wind 'clumping factor' 

The clumping of hot star winds has been the subject of much 
recent interest, as it has become apparent that a clumped wind 
can lead to increased emission for the same amount of mate- 
rial. This means that mass-loss rates derived from spectroscopic 
analysis under the assumption of a smooth outflo w are system- 
atically over-estimated by factors of 3 and higher dBouret et alj 
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120051: iFullerton et all 120061: iPuls et al.|[2006h . To quantify this, 
the clumping of a star's wind is usually parameterized by defin- 
ing the wind clumping factor as a function of the time-averaged 
wind density, 



in the limit where there is negligible inter-clump medium. The 
clumping-factor itself is usually assumed to be a function of 
wind-speed (and hence distance from star), 

/ci = U + (1 - /oo) exp(-v/v d ) (25) 

such that the wind starts off smooth, and clumping is 'switched 
on' at some scale-height above the photosphere when the ve- 
locity reaches v c i. When incorporating clumping into ioniza- 
tion models, some studies have chosen v c i such that clump- 
ing becomes impo rtant at large distances from the star (e.g. 
Dess art et al.l l2000). where hydrodynamical models have shown 
that clumping natu rally arises from instabilit ies in the line- 
driving mechanism dRunacres & Owoc ki 2002). More contem- 
porary studies have chosen smaller values of v c i to initiate 
clumping just above the photosphere, based on the effect a give n 
clumping-law has on spectral line p rofiles (Hil lier et al.l 12003). 
More recently still, Puis et al.] (120061) presented an investigation 
into the clumping factor of different regions of the wind us- 
ing simultaneous optical, IR and radio observations. They find 
that, for stars with winds dense enough to produce line-emission, 
the wind is strongly clumped in the inner 2R+, where the line- 
emission is produced. 

A direct determination of the clumping-law is not possible 
from polarimetry alone. The polarization produced by a clump 
depends on the amount of light the clump intercepts, and the 
thickness of the clump is a redundant parameter provided the 
column-density remains constant. Hence, polarimetry cannot 
distinguish between a dense, flattened clump and a less-dense, 
extended clump, as both have the same optical depth. 

However, the aspects of wind-structure probed by polarime- 
try are useful in conjuction with other observations in determin- 
ing the filling-factor of the wind. It has been shown here that 
polarimetry is sensitive to the number of clumps that the wind is 
broken up into. This may be relevant to other recent studies of 
wind-structure, in particular the study of the X-ray line-profiles. 

Luminous, O-type superg iants have been observed to emit 
X-ray line-emission (e.g. Cas sinelli & Swankl 1 1 983b . likely 
due to shocked gas embedded in thei r outflowing winds 
dLucv& White! [I98(it lOwocki et al.l[l988h . In models of such 
features the red-shifted emission from the back-hemisphere, 
which has a larger path-length to travel to the observer, 
is predicted to be attenuated by t he intervening material , 
producing asymmetric line profiles (lOwocki & Cohenl [2001; 
llgnace & Gavlevll2002l) . However, this is not observed, with the 
X-ra y line profiles of Osgs remarkably symmetric in appearance 
(e.g. lSchulz etaD200(AlWaldron & Cassinellill2001l) . 

In order to explain this phenomenon, it has been suggested 
that either the mass -loss rate must be reduced by a factor of ^ 5 
dCohen et al.l [2006). or that the wind must be clumped in such 
a way as to make the overall wind very porous, i.e. there must 
be large s eparation between the clumps for the photons to leak 
through (lOwocki & Coh en 2007). The porosity-length, h, is de- 
fined as the ratio of the clump size to the volume filling-factor, 
so having a large value of h sets the inter-clump distance to be 
very large compared to the size of each clump. In this regime, 



the X-rays emitted from the back hemisphere can travel through 
the wind relatively unattenuated, and a symmetric line-profile is 
recovered. 

Polarimetry can be used to constrain the available parameter- 
space of these solutio ns. In the porous-wind model of 
lOwocki & Cohenl d2007l) . there is degeneracy between clump 
size and filling-factor - that is, the wind can consist of a few, 
large clumps with a large volume filling-factor; or a large num- 
ber of small, dense clumps and very small filling-factor. These 
two regimes would produce very different polarimetric signa- 
tures. In a wind containing very few clumps, the polarimetric 
evolution of individual clumps could be traced, providing there 
were only a handful of clumps in the inner ~ 5R*. Resolving the 
dynamical evolution of individual clumps would imply that the 
number of clumps was small, and we could then place a lower- 
limit on the filling-factor for a given porosity-length. 

If significant intrinsic polarization was detected, but no in- 
dividual clumps were resolved, this would allow us to place a 
lower-limit to the number of clumps, and hence an upper-limit to 
/. Were no intrinsic polarization to be observed at all, this would 
imply from our models that the number of clumps so large that 
the wind was tending towards a smooth outflow, which would 
not be consistent with asymmetric X-ray line profiles. Given that 
(P) scales with M in our models, such a situation may instead be 
consistent with a lower mass-loss rate for that object. 

6. Summary 

We present an analytic investigation into the wind clumping of 
hot stars, from the point-of-view of their polarimetric variability. 
In particular, we study the effect of varying stellar parameters, 
with the aim of reproducing the level and variability observed in 
LBVs. 

The parameter space was defined such that the only key in- 
put parameter for a given set of stellar/wind parameters was 
the ejection timescale. It was found that, to quantitatively re- 
produce the polarimetric variability of hot stars, the ejection 
timescale must either be very short (<k 1 per hour) or very long 
(~ few per year). In the latter case, the polarization arises from 
the sporadic ejection of very dense clumps, and the variability 
timescale is very long. In the former, the inner wind consists 
of thousands of low-density clumps, with the observed polar- 
ization resulting from random statistical deviations from spher- 
ical symmetry in a slightly-fragmented wind. Here, significant 
changes are expected in night-to-night monitoring. As the po- 
larization of hot stars is observed to be variable on very short 
timescales, and as they are not associated with any significant 
spectroscopic variability that would arise from dense clumps, the 
slightly-fragmented wind explanation is favoured. 

As polarization scales linearly with mass-loss rate, it is spec- 
ulated that all hot stars with high mass-loss rates are likely to 
displ ay polarimetric va riations. This is consistent with the results 
of dDavies et al.ll2005T) which reveal that intrinsic polarization is 
more likely to be found in those LBVs with the strongest Ha 
emission, which is a rough indicator of mass-loss rate. It is also 
consistent with the lower incidence of polarization in O super- 
giants and WRs, which likely have mass-loss rates factors of ^5 
lower. 
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